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We utilize a recent formulation of a spherically symmetric spacetime endowed with a general 
decomposition of the energy momentum tensor [Phys. Rev. D 75 024031 (2007)] to derive equa- 
tions governing spherically symmetric distributions of electromagnetic matter. We show the system 
reduces to the Reissner-Nordstrom spacetime in general, spherically symmetric coordinates in the 
vacuum limit. Furthermore, we show reduction to the charged Vaidya spacetime in non-null coor- 
dinates when certain equations of states are chosen. A model of gravitational collapse is discussed 
whereby a charged fluid resides within a boundary of finite radial extent on the initial hypersurface, 
and is allowed to radiate charged particles. Our formalism allows for the discussion of all regions 
in this model without the need for complicated matching schemes at the interfaces between succes- 
sive regions. As further examples we consider the collapse of a thin shell of charged matter onto a 
Reissner-Nordstrom black hole. Finally, we reduce the entire system of equations to the static case 
such that we have the equations for hydrostatic equilibrium of a charged fluid. 

PACS numbers: 04.40.Nr, 04.20.Ex, 04.20.Jb 



I. INTRODUCTION 

Recently Lasky and Lun [l[ reduced Einstein's field 
equations governing a spherically symmetric distribution 
of a completely general decomposition of the energy- 
momentum tensor to a system conducive to an initial 
value problem. It was shown that under certain equa- 
tions of state, the system reduced to the Vaidya space- 
time in non-null coordinates. Furthermore, when all 
energy-momentum terms vanished the system reduced 
to the Schwarzschild spacetime in general coordinates. 
This formulation enables all regions of the spacetime to 
be written in a single coordinate system with equations of 
state tending to certain values at the interfaces between 
the regions. Therefore, the use of junction conditions at 
the interfaces is not required as one simply establishes 
an initial hypersurface and lets the system evolve via the 
governing set of equations. 

The presentpaper highlights the generality of the fluid 
elucidated in [1J by incorporating electromagnetic (EM) 
fields into the spacetime. In this way we are doing plasma 
physics in full general relativity and we show that the 
equations have familiar Newtonian analogues. The role 
of EM fields in fully relativistic gravitational collapse has 
been the topic of much interest over recent years @, H, H| • 
Of particular interest is the formation of naked singular- 
ities [5|, [6|, [7[, the formation and evolution of charged 
neutron stars within the context of full general relativity 
[1, [t| and the collapse of charged shells of matter 
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onto existing black holes [11] . We utilize a formalism 
which puts general, spherically symmetric charged fluids 
into an initial value setting on a spacelike hypersurface 
such that all of these issues can be discussed either ana- 
lytically or numerically. This paper concentrates on the 
formulation of the plasma equations, with an emphasis 
on the reduction of the full set of equations to various so- 
lutions of charged fluids in spherical symmetry, including 
the investigation of the collapse of infinitely thin shells. 

Known charged fluid solutions in general relativity in- 
clude the charged Vaidya and Reissner-Nordstrom (RN) 
spacetimes. The generalized BirkhofFs theorem implies 
the RN spacetime is the unique, spherically symmetric 
vacuum solution endowed with an electric field [l2j . The 
charged Vaidya solution [ID, 0, HH, [H| is a generaliza- 
tion of the Vaidya metric which describes null radiation. 
Both the Vaidya and charged Vaidya spacetimes are usu- 
ally exhibited by using null coordinates, and indeed the 
mathematical form of the metric in these coordinates is 
rather simple [UGH- However, the physical interpreta- 
tion of the fluid and the matching of the fluid to other 
spacetimes is fraught with difficulties [HI, |20| . In [l[, 
we re-cast the Vaidya spacetime in non-null coordinates 
such that it is conducive to an initial value formulation 
and showed that this made the matching of spacetimes 
significantly simpler. In this article we extend the result 
to include the charged Vaidya spacetime in spacelike co- 
ordinates and further discuss the matching of this to an 
interior charged fluid and also to an exterior RN solution. 

Another example we elucidate in this paper is that of 
an infinitely thin charged shell falling onto an RN black 
hole. We are able to describe this entire scenario in terms 
of a single line element and hence as a single solution of 
the Einstein field equations. Furthermore, the differen- 
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tial equation governing the collapsing charged shell is re- 
duced to quadrature, implying the dynamics of the entire 
spacetime is known. 

The Einstein-Maxwell equations describe the coupling 
of a gravitational field to an EM field. These comprise the 
Einstein field equations whereby the energy momentum 
tensor has a contribution from the Faraday tensor. The 
resultant energy momentum tensor can be represented as 
simply the sum of two fields [T^j 



satisfies Frobenius theorem, 



n[^V v n a ] = 0, 



(2) 



where square brackets denote antisymmctrization and 
V M is the unique four-metric connection, implies the nor- 
mal vector is hypersurface forming. Further define a pro- 
jection tensor, _L Mt/ , which projects all quantities onto the 
spacelike hypersurfaces formed by the normal vector, 
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Throughout the article a superscript "F" represents en- 
tities relating to the fluid contribution to the energy mo- 
mentum tensor, T^ v , and a superscript "EM" relates 
quantities to the EM contribution to the energy momen- 
tum tensor, T^ V M . In section UU of this article we review 
the form of T^ V M and the Maxwell equations in a co- 
variant formulation. Furthermore, in this section we de- 
compose Maxwell's equations and T® V M into spatial and 
temporal components by exploiting the 3 + 1 formalism 
for general relativity. Section IIIII then reviews the neces- 
sary components of the general decomposition of and 
the equations derived in We point out the existence 
of the generalization of the Lorentz force law when EM is 
coupled to the Einstein field equations. This equation is 
often neglected when studying Einstein-Maxwell fields, 
however it is a necessary component of the system of 
equations when modelling physically realistic scenario's. 
Section HVl then puts together the results of the two pre- 
ceding sections to derive the reduced field equations gov- 
erning a spherically symmetric distribution of charged 
matter. In section |V] we show how the system reduces 
to the RN spacetime in the appropriate limit and section 
I VII shows the derivation of the relevant equations of state 
and the reduction of the system to the charged Vaidya 
spacetime. We move to a slightly different application in 
section I VIII to discuss in detail a shell of charged matter 
collapsing onto a black hole. We show the entire system 
is given by a single line element, where the differential 
equation of the shell reduces to quadrature. Finally, in 
section IVIIII we reduce the system to the static sub-case 
and hence derive the equations of hydrostatic equilibrium 
for a charged fluid ball. 

Throughout the article we use a +2 signature and co- 
ordinates are denoted by a; M = (i, r, 6, <p) unless otherwise 
stated. Geometrized units are used, whereby c = G = 1 . 
Greek indices range from ... 3, Latin from 1 ... 3 and 
all other conventions follow 11211. 



II. COVARIANT ELECTROMAGNETISM 

We provide here a brief introduction to the relevant 
components of the 3 + 1 formalism for general relativity. 
Begin by defining a vector field, n^, which is timclike 
and normalized, n a n a = —1. Further demanding this 
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where is the four- metric. 

The energy momentum tensor for EM fields is ex- 
pressed in terms of the anti-symmetric Faraday tensor 
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An observer moving in the direction of the normal vector 
will then measure the electric and magnetic field intensi- 
ties respectively as 



E„ 



F n a 

A fj,a 11 



and 



B„ 
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where e^ a is the three-Levi Civita pseudo-tensor. One 
can show from the anti-symmetric nature of the Faraday 
tensor and the spatial nature of the three-Levi Civita 
pseudo tensor that E a n a = B a n a = 0, implying the 
electric and magnetic fields are spatial. The Faraday 
tensor can be uniquely decomposed onto and orthogo- 
nal to the spacelike hypersurfaces specified by its normal 
vector. After some algebra we can show 
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Using the above information, equation ^ can be irre- 
ducibly decomposed resulting in the energy momentum 
tensor being expressed in terms of the electric and mag- 
netic fields 
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We note the vector jjf M is the EM Poynting vector [22 
which is a measure of the energy flow in the spacetime. 
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Maxwell's equations are expressed covariantly in terms 
of the Faraday tensor according to [2l| 

V [<r F H =0, (12) 
V Q F MQ =4ttJ m . (13) 

Here, is the charge-current four-vector, which can be 
decomposed according to 

i7 M = V + erV (14) 

where := _Lp Q j7a is the three-current density and 
e := —n a J a is the charge density. By taking the diver- 
gence of both sides of equation (fTB"]) , we obtain the conser- 
vation law for the charge-current four-vector, V a J a = 0. 
Decomposing this with respect to the normal implies 

C n e -Ke + {D a + h a ) i a = 0, (15) 

where C n is the Lie derivative operator with respect to 
the normal vector field, K is the trace of the extrinsic 
curvature, Di is the unique three-metric connection asso- 
ciated with _Ljj and :— n Q V Q n M is the observers four- 
acceleration. Equation (|15[) is an identity if Maxwell's 
equations are satisfied. However, it can be used if one is 
given either an arbitrary charge distribution or the three- 
current density to determine the other unknown quantity. 

Decomposing equation ([12")) with respect to the normal 
vector provides two equations 

D a B a =0, (16) 
C n B, t - KB^ + 2K«B a = - e Ma/3 (D a + n a ) E 13 , (17) 

where K^ v is the extrinsic curvature. Decomposing equa- 
tion (TIT?]) gives 

D a E a =4ne, (18) 

C n E^ - KE^ + 2K„ a E a =e^ p (D a + n a ) B 13 - 4ttv 

(19) 

Equations (|16M19[) are the more familiar version of 
Maxwell's equations in curved spacetimc. 

To complete the description of EM in general relativ- 
ity we now derive the Lorentz force law. We expand the 
divergence of the EM component of the energy momen- 
tum tensor, V Q T^ M , with respect to the definition in 
terms of the Faraday tensor Q. Then by using both 
the anti-symmetric nature of the Faraday tensor and the 
four-dimensional Maxwell equations, (fT2"| and (fT3"|) , we 
can show 

VI™ = -F m J a . (20) 

Conservation of energy momentum in general relativity 
is a statement that the divergence of the total energy mo- 
mentum tensor vanishes. By expanding this in terms of 
the EM and fluid components, equation JT]), and substi- 
tuting equation (|20p we find the conservation of energy 
momentum in general relativity gives 

v% F Q = i^r. (2i) 



This is a four-dimensional equation which can be decom- 
posed into three constraint equations on the hypersur- 
face and a single propagation equation. This decom- 
position involves taking the L^V^T^ component and 

the n a V^T^ component respectively. However, the left 
hand sides of these equations requires the explicit knowl- 
edge of the decomposition of the fluid component of the 
energy momentum equation which shall be discussed in 
the next section. 

Summarily, equations (|15H19p provide the full set of 
Maxwell's equations governing EM fields in curved space- 
time. These couple back to gravitation through the con- 
servation of energy momentum tensor given by equation 

(ED- 



III. FLUID FORMULATION 

A. Covariant equations 

In [l[ we used a generic decomposition of the total 
energy momentum tensor given by 

T MV =pn ll n v + P ± MV +2j( /1 n I /j + U^, (22) 

where p and P are the total energy density and total 
isotropic pressure respectively, j M is interpreted loosely 
as a heat flux vector and n Mt/ as an anisotropic stress 
tensor. 

We use the same decomposition for the fluid compo- 
nent of the energy momentum, T F , as we used for the 
total energy momentum tensor, i.e. equation (|22p . How- 
ever, we denote each quantity decomposed from the fluid 
contribution to the total energy momentum tensor with 
a superscript "F"' . The total energy momentum compo- 
nents are now given in terms of their constituent fluid 
and EM contributions by using equations |T]), (|51 fTT|) and 
(1221) . and we find 

p = p F + ±- (E a E a + B a B a ) , (23) 
P = P F + -L (E a E a + B a B a ) , (24) 

/47T 

j^=j^ + ^s fiaP E a B' 3 ) (25) 
H^=H^+n^ M , (26) 

where n^j 11 is given by equation (fTTj) . 

Before moving on to the spherically symmetric reduc- 
tion of the equations, we discuss the conservation of en- 
ergy momentum equations which lead to a generalized 
Lorentz force law. As mentioned above, we decompose 
equation (|20[) onto and orthogonal to the spacelike hyper- 
surface formed by the normal vector. The first equation 
can be found from the -L^V^T^g part of the equation, 
which after some algebra and using equations ©, (fT4"l) 
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and (|22p can be shown to give 



+ (P F 



n a )Il F a +D^P F 



P* 



= eEfj, 



a B f) . (27) 



This is a covariant representation of the Lorentz force 
law coupled to the gravitational field, which can be seen 
intuitively by considering an example of the flow of some 
current. Here, describes the flow of the charge and 
j F describes the flow of the "fluid particles" . As the 
charge is tied to the fluid, these two quantities essentially 
describe the same phenomena. Therefore, the term C n j F 
describes a timelike derivative of a flow of fluid particles, 
i.e. the rate of change of momentum of fluid particles. 
Furthermore, the right hand side is the electric field plus 
the wedge product of the three-current density and the 
magnetic field. Feynman [23| states that this force law 
is required in classical EM theory because having all the 
electric and magnetic fields doesn't tell us anything until 
we know what they do to the charges (see chapter 18-3 of 

1H). 

The second term from the decomposition of the conser- 
vation of energy momentum is given by n a V^T ( ^g. Again 
after some manipulation this becomes 



r I 



P F ) K+ (D c 



)3a 

-a/3 _ 



n:^ Q " = -E a i c 



(28) 



The left hand side of this equation is equivalent to the 
continuity equation of fluid dynamics. The term on the 
right hand side describes the electrical energy for the 
plasma [2~4j . This equation has a Newtonian analogue in 
plasma physics as an alternative form of the heat equa- 
tion and it represents the fact that the gain in energy of 
the material is due to many different facets including heat 
flux, electrical energy and the work done by the pressure 
(see equation 2.40 of Priest HH). 



B. Spherical Symmetry 

Spherical symmetry implies the density and pressure 
are just functions of the temporal and radial coordinates. 
We can also write the heat flux vector as 

J„ = W(t,r), j(t,r),0, 0]. (29) 

Furthermore, we define a spatial tensor according to 

P/ :=diag[-2, 1, 1]. (30) 

Now, the anisotropic stress tensor, is symmetric, 

spatial (ILj Q n Q — 0) and trace-free (II Q a = 0). It is 
trivial to show with spherical symmetry we can always 
write 



rV : = II(t,r)iV, 
where H(t, r) is the distinct eigenvalue of Hy. 



(31) 



In [l[ we further began with a spherically symmetric 
line element expressed in 3 + 1 form 



dS 2 = -a 2 dt 2 



(f3dt + dr)' 
1+E 



R 2 dfl 2 , (32) 



where a(t, r) > is the lapse function, /3{t, r) is the ra- 
dial component of the shift vector, E(t,r) > —1 reduces 
to the Lemaitre-Tolman-Bondi energy function in the un- 
charged dust limit [25[ and dfl 2 := dd 2 + sin 2 6d(j) 2 . Fur- 
thermore, R = R{t 1 r) is an arbitrary function whereby, 
without loss of generality, we set R(t, 0) = 0. The normal 
vector in these coordinates was shown to be given by 



n" = -(-l, (3, 0, 0). 
a 



(33) 



We were then able to show the equations governing the 
system were given by the following complicated set of 
equations. Firstly, we defined a "mass" function in terms 
of the total energy density, p, and the total heat flux, j, 
by 



dM ( dR .„ , , r , 

— :=4irlp— +j£ n R)R 2 



(34) 



Equations governing the evolution of E, R and M were 
then found to be given by 



dR C n E ldRdfS d f 

to WT^) = 4njR -a^d^-Yr {£nR} ' (35) 



(«) 2 ^-l + (l + *)(f 



-AttR 2 



(P-2U) C n R + j^(l + E) 
or 



(36) 
, (37) 



M , , 1 + EdadR „ 2 

— + 4tt P - 211 )R = — — - £„ 2 P, 

R z a or or 



(38) 



where C n R := C n (C n R). Furthermore, the final equa- 
tion (|38p is often trivially satisfied by the three preceding 
equations, however there are circumstances where this 
equation is required as we will see below. There are two 
further equations coming from the conservation of energy 
momentum, equations (|2"T|) and (j2"B"|) . Using the metric 
(|32)l to express these equations they become long and aes- 
thetically awkward and for this reason we suppress the 
independence of the EM and fluid contributions for the 
energy momentum tensor. These two equations become 
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respectively 



gn (j\/T+£) 



2j 



1 13/9 
2(l + £) " ad^ 



i? 



( P + p-2U) da d | enaa 

a 9r 9r i? dr 



(39) 



£„p = (p + P - 211) 



1 



j(l + E) 



ldj_ 

3 dr 



2(1 + E) 

1 dE 
2(l + E)~dr~ 



15/3 2 

a S>-nR 

a or R 



2 8R 
R~dr~ 



2 da 
a dr 



(40) 



The system requires two further equations of state, for 
example where P and II are expressed in terms of p and j . 
We will show in the coming sections that equations (13"9")) 
and (|40[) will reduce considerably for various subsystems 
of the general equations. The two equations above are 
equivalent to the ones discussed at the end of section HT1 
If we expand these equations in terms of the fluid compo- 
nents and the EM components of the energy momentum 
tensor, then indeed the term on the right hand side of 
equation (j2"?| appears in equation (f3"9"]l . 

Summarily, the entire system of equations governing 
a spherically symmetric distribution of fluid is given by 
the line element (32) . where the metric coefficients are re- 
lated to the total contribution of the energy momentum 
tensor through equations (|34ll40p . This total energy mo- 
mentum tensor is decomposed into constituent fluid and 
EM contributions by equations (|23H26[) . In the next sec- 
tion we derive the electric and magnetic field components 
in spherical symmetry using the metric (|32[) . 



IV. CHARGED FLUID 



spacetime, which would intuitively introduce a non-zero 
magnetic field. However, as the current in the space- 
time must be uniform across the two-spheres the mag- 
netic field created everywhere will exactly cancel itself 
out. Expressing equation (|17p in terms of the metric co- 
efficients we find it vanishes identically. Equations (fTS)) 
and (|19p respectively become 



ER 2 



0_ 

dr VVl + E 
£R 2 



r 

*~"n 



4neR 2 
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-AiriR 2 



VTTe J VTTe' 



(46) 
(47) 



Integrating equation (|46|) we find the electric field is given 
by 



c _ QVTTE 
*~ R 2 ' 

where we have defined Q(t, r) such that 

f r eR 2 
Q := 47r / dr 
Jo VTTe 



(48) 



(49) 



This is interpreted as the total charge within a two-sphere 
of radius r 2}. We will see in the vacuum Einstein- 
Maxwell case 1 , it reduces to the charge of the RN so- 
lution (see section fV]) . Putting (|4"8")) back through the 
remaining Maxwell equation (|47|) . we find 



CnQ = 



-4iriR 2 

VTTe' 



(50) 



Equations (|45|) and (|48|) give us equations of state for 
the EM contribution to the energy momentum tensor. 
This is seen by evaluating equations (|51fl"Tj) with respect 
to Q, from which it can be shown 



Spherical symmetry implies all vector fields can be 
written with only non-zero radial components. There- 
fore, we make the following definitions 



E"=[0,£(t,r),0, 0], 
B M = [0, B(t,r), 0, 0] , 
i" = [0, i(t,r), 0, 0] , 



(41) 
(42) 
(43) 



where i is not to be confused with the complex number. 
Maxwell's equations can now be expressed using the line 
element (|3"2")l . Beginning with equation (|16[) . we find 



d ( BR 2 



dr \ VTTE 



=0. 



(44) 



Integrating this and using R(t, 0) = implies the mag- 
netic field in our frame vanishes everywhere 



EM 



EM 



J 



EM 



J R EM 

2 
=0. 



Q 2 



87T.R 4 ' 



(51) 
(52) 



As mentioned, the j EM component of the Maxwell energy 
momentum tensor is the Poynting flux. This term being 
zero implies there is no transference of energy via the EM 
fields which is a direct result of the spherical symmetry. 
As this term is zero, for the remainder of the article we 
drop subsequent superscripts on the j F term. 



Finally, equation (|15p can be expressed in terms of the 
metric coefficients as 



In 



/ eR 2 



\VTTe 
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In 



iR 2 a 

VTTe 



(53) 



B = 0. 



(45) 



This is a direct consequence of the symmetry of the prob- 
lem. We note that we can still have a current in the 



1 Throughout the article we refer to a vacuum in the sense that 
the fluid component of the energy momentum tensor vanishes, 
i.e. = 0. Of course in this case the total energy momentum 
tensor is still non-zero due to the EM component. 
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We have essentially shown that in the spherically sym- 
metric Einstein-Maxwell system, the Maxwell equations 
can be integrated completely. On the initial spacelike hy- 
persurface we stipulate the charge density and the three- 
current density such that equation (|53[) is satisfied. Fur- 
thermore, the Einstein constraint equations can be solved 
such that we know the geometry of the initial spacelike 
hypersurface, implying the metric coefficients are known 
initially. Therefore, through equation (|49|) the charge 
per unit volume, Q, is known on the initial hypersur- 
face. Equation (|50|) can then be integrated to give Q on 
subsequent spacelike hypersurfaces. 

Using (T5"Tj) , ([52")) and the definition for the mass func- 
tion (l34l). we can show 



8M 
dr 



47T 



,dR 
dr 



jC n R)R 



and equation (|57|) becomes 



4irR 2 
Q 2 



Q 2 dR 



OR 



(54) 



(P F -2U F ) C n R + j—(l + E) 



2R 2 



C n R. 



(55) 



The following sections concentrate on reducing the 
general system of equations to specific examples and 
models. 



V. EXTERIOR VACUUM REGION 

A vacuum region is defined such that all fluid compo- 
nents of the energy momentum tensor vanish, i.e. 



p F = P F = f = U F = 0. 



(56) 



Furthermore, if we are in vacuum there are no charge 
carrying particles and hence 



(57) 



s = i = 0. 



Using this we can integrate equations (|49|) and (|50|) to 
find 



Q = Qo, 



(58) 



where Qo is a constant equal to the total charge inside a 
sphere of certain radius being the boundary marking the 
vacuum region. The mass function from equation (|54[) 
can now be integrated and we find 



M = M - 



2R 



(59) 



where Mq is a constant of integration. 



There are three remaining non-trivial equations, (I35p . 
([56]) and which become respectively 



C n E 1 df3 
2{l + E) + a~dr~ 



dr 



(C n R) 2 -(! + £) 



1 + E da OR 
a dr dr 



dR 
dr 



{CnR) , 

2M Q , Qo' 



R 



R 2 



Mo 



Qo 2 



r 1 r — — _ 



(60) 

)• 

(61) 
(62) 



The generalized Birkhoff's theorem states a spherically 
symmetric vacuum with an electric field is necessarily the 
RN spacetime [T^]- Therefore, the metric (|32|) where the 
coefficients are solutions of equations (|60ti62|) must be the 
RN spacetime. The amount of generality corresponds to 
the vast generality for the interior regions of the space- 
time. Upon determining an interior, charged fluid filled 
region, one simply lets the metric coefficients be contin- 
uous across the boundary of the interface to determine 
the specific form of the RN region. 

It is trivial to see the equations governing the RN 
spacetime are a generalization of the equations govern- 
ing the Schwarzschild spacetime given in [1| , whereby the 
Schwarzschild set are retrieved simply by setting Qo = 0. 

In general, the solution of equations (|60H62[) may be 
non-trivial. However, we can pick out a handful of exact 
solutions by making ad hoc stipulations for some of the 
metric coefficients by using residual coordinate freedom 
available. For example, by setting f3 = and R = r, 
equation (|6Tj) gives E in terms of Mo and Qo and equa- 
tion (|6"0"|) can subsequently be integrated for a and the 
metric is expressed in the usual RN coordinates. Alter- 
natively, Novikov-type coordinates are found by setting 
a = 1 and (3 = 0. This implies equations ([60]) and (f6"Tj) 
can be integrated such that the line element is 



dS l 



-dt 2 



1 (dR 



T\ V dr 



dr z + R 2 dn 2 



(63) 



where Ti(r) is an arbitrary function of integration and 
R is given implicitly by the solution of 



dR 



-l + 2Mo/R-Qo 2 /R 2 +Fi{r 



(64) 



where J^2{ r ) is another function of integration. The so- 
lution of this integral can be written down, however it is 
complicated and provides no new insight for the present 
discussion. 

Another set of coordinates can be obtained by assum- 
ing no temporal dependance in the line element. Also 
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letting R = r implies the line element becomes 



dS 2 



2Mn 



Go' 



dt 2 



E + 2M /r- ql/r 2 ,,, 
±2 V 1 + 25 

dr2 2 ,^2 

+ TTE +rdQ - 



(65) 



This line element still contains a single free function of 
the radial coordinate, E(r), which may be specified ar- 
bitrarily. We note it is equivalent to the general line 
element expressed by Patino and Rago [26j, and if one 
chooses (1 + Ey 1 = 1 - 2M /r + Q 2 /r 2 , then without 
loss of generalit y, it results in the line element expressed 
by Papapetrou [27j . 

There is one more coordinate system that will prove 
relevant for section I VIII on the collapse of infinitely 
thin shells and thus we introduce it here. Generalized 
Painleve-Gullstrand type coordinates (2f| have canonical 
two-spheres and zero acceleration, i.e. R = r and a = 1 
respectively. Putting this through equations (|60H62[) wc 
find 



/? 2 



E 



C n E = 0, 
2M 



(66) 
(67) 



Therefore, the only free function in the metric is the en- 
ergy function E, which is a solution of (|66p . Furthermore, 
when taking the square root of (|67[) we can either choose 
a plus or minus sign, which give different regions of the 
RN solution on the Kruskal diagram (see for e.g. [12]). 
A trivial solution of (|66p is E = 0, which are Painleve- 
Gullstrand coordinates for the RN spacetime [H, [29[ (see 
[30l ] for a review of these coordinates in the Schwarzschild 
case) . For an expose of more solutions of these equations 
and their properties see (3lj . 



VI. CHARGED VAIDYA 

In [l[ we used a coordinate transformation to go from 
the outgoing Vaidya metric to the 3 + 1 line element, 
equation [[35]). Here we generalize that result to include 
the charged Vaidya spacetime [HI, [3, [H, [l(| (also see 
[l7l Il8j for more recent results). We begin with the out- 
going charged Vaidya metric 

dS 2 = - (l - ^) dv 2 - 2dvdR + R 2 dVL 2 , (68) 

where v is an outgoing null coordinate, (R, 9, <p) are 
spacelike coordinates and m — m(v,r). The energy mo- 
mentum tensor may be decomposed with respect to two 
null vectors [il| 

-2]uu — P ^fl^is 

+ (w + a) {ifjtu + kpt v ) + crg^u, (69) 



where the Einstein equations can be shown to satisfy 
P 1 dm 



AttR 2 dv ' 
1 dm 

AttR 2 M' 
1 d 2 m 

8ttR dR 2 ■ 



(70) 
(71) 
(72) 



Furthermore, the vectors in (|69p are given in component 
form by 



/ —t) v 



(73) 
(74) 



One can see in (1691) if lu = a = 0, then the energy momen- 
tum tensor reduces to that of null dust and the spacetime 
reduces to Vaidya, whereby p F is the only non-zero com- 
ponent of the energy momentum tensor. As p F in the 
Vaidya system is the only contribution to the energy mo- 
mentum tensor which describes an uncharged null fluid, 
it can be shown lu and a contribute to only the EM part 
of the energy momentum tensor in the charged Vaidya 
system [3. 

We perform a coordinate transformation from the 
charged Vaidya line element to the non-null metric (|3"2"|) 
such that v = v(t,R) and R = R(t,r). This coordinate 
transformation is given by 



OR am (3(1 -m/R) OR -(l-m/R) 

~dt ~ ~~R VTTe ' ~ VTTe 

dv (3 dv 1 

~dl ~ a+ VT"kE' dr' ~ VTTW' 



(75) 



We also write here the inverse coordinate transformation 
as this will be necessary for the later derivations 



dt -1 / m \ dt -1 
dv a V RJ dR a 



dr 8 / m\ m , — Or 13 /- 

dv a \ RJ R dR a 



dr 13 



E. 
(77) 



One can see the coordinate transformation (|75[) implies 
the relation 



C„R 



rn 
~~R~' 



(78) 



Putting this through equation (|36|) we find the result that 
the mass defined in equation (f3"4"]l reduces to the mass of 
the charged Vaidya system 



M = m. 



(79) 



Now, for solutions of the differential equations defining 
the coordinate transformation to exist and hence for the 
coordinate transformation to be valid, the integrability 
conditions of the coordinate transformation must be sat- 
isfied. These integrability conditions provide equations 
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constraining the form of the metric coefficients, 
some algebra we find these to be 



VT+Edc 



C n E 



dr 2(1 + E) 



180 

a dr 

f VTTe 



After 



(80) 



~dr~) 



M 

-w (81) 



Putting the integrability conditions through equations 
(|35[) and (|38|) we find after some algebra 



£ n M+VTT~E 



8M 
dr 



= 47rir 



p f -2n F -jy/T+E)-^. 



(82) 



Furthermore, putting the integrability equations and the 
coordinate transformation through equations (|34[) and 
(|37|) and substituting through equation (f82|) , we can show 
the fluid variables satisfy a general equation of state 



P l 



21T 



2 j VTTe = o 



(83) 



Now, by expanding the derivative with respect to the 
null coordinate v in equation (|T0[) with respect to the 
temporal and radial coordinates, t and r, one can show 



dM 
dv 



= -47ri? 2 jVf + E. 



(84) 



Hence, direct comparison with equation (|70[) implies 
" E = .iVl 



P 



■■ jVl + E 



(85) 



Doing the same for equation (|TTj) and (f?2")) we find re- 
spectively 



Q 2 



87T.R 4 

2R 



e — 



Q 2 



VTTe) 8vri? 4 ' 



(86) 
(87) 



Now, substituting the conditions for the coordinate 
transformations and equations (|83H87[) into the conserva- 
tion of energy momentum equations, (|3"9"j) and (|4D)) . after 
much algebra we can show these reduce to the simple 
relation 



P EM = P F 



31F 



The uncharged Vaidya system had four energy momen- 
tum variables with three equations of state implying we 
could express all energy momentum variables in terms of 
just one of the variables. However, the charged Vaidya 
system has six energy momentum variables, p F , P F , H F , 
j, p EM and a. However, we still only have three equa- 
tions of state which can be summarized from the above 



as 



=jVT 
=p F - 



-E 
2U F 



-p EM + a 



3W 



(89) 
(90) 
(91) 



It is trivial to show when p EM = a = 0, the system re- 
duces to the uncharged Vaidya in non-null coordinates 
The extra generality in this system is attributed to 
the extra generality associated with the charge distribu- 
tions and currents throughout the spacctime. 

Two fmore equations from this system will prove useful 
in the further analysis. Firstly, putting the equations of 
state through equation (|52")) we find 



C n M + VTTe^- 

or 



Q 2 
2B? 



iw>- (92) 



Secondly, putting this equation back through the inte- 
grability condition, equation (f8Tj) . implies 



I da 

a dr 



-I 



R 2 VT 



E \ 2R 



(93) 



Summarily, the system of equations governing the 
charged Vaidya spacetime in non-null coordinates is the 
line element (|32[) where the metric coefficients are given 
by dill), (HO]), (|92]1. and {93]). One can imagine establish- 
ing an initial value problem such that an initial spacelike 
hypersurface is established with an interior charged fluid 
ball reaching some finite radius and a vacuum RN region 
extending beyond this radius. The system then evolves 
and the fluid ball is allowed to radiate null radiation such 
that a region of charged Vaidya as described above devel- 
ops between the fluid ball and the vacuum. Taking any 
subsequent t equals constant slice of the spacetime will re- 
sult in a region of a charged fluid ball with charged Vaidya 
extending beyond this and RN reaching out to radial in- 
finity. In this way causality is provided for the charged 
Vaidya spacetime within a realistic physical model. Of 
course the equations described herein are difficult to solve 
analytically, as such future work will provide a numerical 
realization of this paradigm. 



VII. THIN CHARGED SHELL 

Consider an infinitely thin shell of charged matter 
falling into an RN black hole. As the shell is infinitely 
thin, there are no internal forces between particles other 
than the Coulomb forces implying the fluid is dust, i.e. 
P F = H F = j = 0. Therefore, the only non-zero contri- 
bution to the fluid component of the energy momentum 
tensor is from the energy density, which is given by 



p F S(r)+p F S(r-n) 



(94) 



Here, S is the Dirac delta, p F € R is the energy density 
of the initial RN black hole situated at r — and p F € K 
is the energy density of the shell, which is located at r = 
r*(t). Furthermore, the dust particles of the collapsing 
shell carry charge and hence the charge density is given 
by 



e 6(r) + £+8 (r - r*) , 



(95) 
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where £q g M is the charge density of the point mass at 
r = and <E K is the charge density of the shell. 

The energy momentum tensor here is described by two 
variables, p F and e. As there are four degrees of freedom 
in the metric (a, 3, E and R) and only two degrees of 
freedom in the energy momentum tensor, we have resid- 
ual coordinate freedom in the line element (see the dis- 
cussion in 1]). Hence, without loss of generality we use 
canonical two-spheres in the following derivations, i.e. 



R = r. 



(96) 



We will use the remaining coordinate freedom associ- 
ated with the spacetime implicitly at a later stage of the 
derivation to further simplify the mathematics. Equation 
(|4"9")l can now be integrated to give 



Q= Qo + QM (r - r*) 



(97) 



where Qq, Q* € K are constants of integration which rep- 
resent the charge at r = and of the shell respectively 
and U is the Heaviside step function. Integrating the 
definition of the mass function (l34l) we find 



M =M - 



2r 



M* - (2Q + &,) 
2r 



M -% 



2r 



r < r* 



(M + M*) 



(Qa + Q.f _ > 
2r 



U{r- r*) 
(98) 



Here, Mo is the mass of the initial Schwarzschild black 
hole and M* is another constant of integration which is 
the mass of the shell. Equation (|9"8")) gives the expected 
result that for r < r* the spacetime is the RN spacetime 
with mass and charge contributions being given by M 
and Qo respectively. The region for which r > r± is also 
the RN spacetime, however the mass and charge contri- 
butions are given by Mo + M* and Qo + Q* respectively. 
This further implies that if the charge at r = and at the 
shell are of equal but opposite sign (for example protons 
and electrons), then the exterior region of the spacetime, 
r > r*, is given simply as the Schwarzschild spacetime 
with total mass Mo + M*. That is, exterior to the shell 
in this scenario one can not feel the effects of the charges 
(this is discussed more in section [VII C 3|) . 

Now, putting the derived form of the mass function 
through equation (|37[) . we find after some algebra 



M 



0* 

2r 



(2Q + Q*) 



C n [U (r - r*)] = 0. (99) 



Hence, provided M ^ (2Q + Q*) /2r, we find 

dr* 
dt 



+ /3U(r-r*) = 0. 



(100) 



This is trivially satisfied for all r r+ and provides no 
new information. However, for r = r t , we see 



dt ~ P *' 



(101) 



where /3* := 3 [t, r*(t)] is the shift function evaluated on 
the shell of matter. Equation QlOip implies the observer 
is moving with the velocity of the shell of matter. 

Evaluating the two conservation of energy momentum 
equations, we find (|40|) is trivially satisfied and ([39]) be- 
comes after some algebra 



Po 5 'tf(r) + pi '8 (r - r*) da Q* 



a 



dr 87rr 4 



{2Q + Q*)6{r 



-n). 

(102) 



This equation can be integrated directly and one recov- 
ers two boundary conditions for the lapse function. The 
first of these boundary conditions is at r = 0, however 
there is a physical singularity at this point implying the 
spacetime is only well defined for r > 0. Therefore this 
first boundary condition is physically irrelevant for the 
system in study. The second boundary condition derived 
from the integration of equation (|102|) is at the boundary 
of the shell of charged matter 



where we have defined the constant 



k := 



g* (2Q + g*) 



(103) 



(104) 



This boundary condition implies the lapse is cither 
steeply increasing or decreasing across the boundary of 
the charged shell depending on the sign of the respective 
charges. Furthermore, we see that if the charge of the 
shell, Q+, and the charge of the point mass, Qo, are of 
the same sign then the gradient of the lapse at is nec- 
essarily increasing. The case where k = is a greatly 
simplified scenario as it implies a can be set to unity 
throughout the spacetime (see section IVIIC 21) . 

The lapse and energy functions are related through 
equation (JTSJ) , which reduces to 



2{l + E) 



P d 

-— (In a) . 

a or 



(105) 



This equation can be evaluated on the shell and we find 
by using equation (jlOip and the chain rule that the left 
hand side simply becomes an ordinary timelike derivative 
of the energy function along the shell. Furthermore, by 
using the boundary condition (|103[) . the entire equation 
can be re-expressed as 



d_ 

dt 



(W 



l + £* 



k dr+ 
ri~dt' 



(106) 



where E+(t) := E [t,r*(t)]. Upon integrating this equa- 
tion with respect to t we introduce an arbitrary constant 
which relates the initial energy function on the shell to 
the initial radius of the shell. We can arbitrarily stipu- 
late the initial energy function on the shell such that this 
arbitrary constant equals unity without loss of generality. 
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It is subsequently found that the energy function on the 
shell is given by 



1 + .E* 



exp 



-2k 
3r2" 



(107) 



As mentioned, E+ is the value of the energy func- 
tion along the shell as it evolves through the space- 
time. Therefore, in keeping consistent with the current 
paradigm, we choose to let E(t, r) be a continuous func- 
tion throughout the spacetime such that for r = r+ the 
energy function reduces to the form given in equation 
(I107p . That is we stipulate E(t,r) for the entire space- 
time to be 



1 + E = exp 



-2k 
"3^3" 



(108) 



We note here that this actually implies the energy func- 
tion is independent of the temporal coordinate. In doing 
this we have lost some generality in the solution as to keep 
the solution as general as possible would have required 
the introduction of extra arbitrary functions whilst still 
ensuring the form of E reduces to that of E+ at r = r*. 
However, the choice we have made above allows the solu- 
tion of the lapse function to be solved with ease. We will 
show that this loss in generality is equivalent in the non- 
charged case to choosing the marginally bound solution, 
i.e. E = 2 . Furthermore, this is equivalent to choosing 
the initial velocity of the shell on the initial hypersurfacc. 

By putting the form of the energy function, equation 
(llOSp . back through equation (| 1 05[) we find the lapse 
function is also independent of the temporal coordinate. 
This simplifies the mathematics as it implies the shift 
function can be cancelled and the equation can be di- 
rectly integrated for the lapse function, which we find 
becomes 



a 



exp 



-k 
3^3 



(109) 



Here, an arbitrary constant of integration has been set to 
unity by re-scaling the temporal coordinate. It is trivial 
to show that this form of the lapse function satisfies the 
boundary condition given by equation (|103D . 

Finally, equation (|3l)|) with R — r implies the shift 
function can be expressed in terms of the lapse, energy 
and mass functions according to 



(3 = a \E 



(110) 



Taking the positive square root of this equation gives a 
model which is collapsing and the negative root gives 



2 The marginally bound solution does not exist in the charged case 
as if we let E = in equation l|105j l then in general the boundary 



an expanding shell [25(. We now know all the metric 
coefficients for the spacetime and hence the system is al- 
most entirely determined. The only remaining piece of 
information is to determine the position of the shell as a 
function of the temporal coordinate. We note that equa- 
tion (|110p can be evaluated on the shell and this gives 
a first order, non-linear ordinary differential equation for 
the position of the shell as a function of time 



dr+ 
dt 



= - e -fc/ 3 ^ 4 e -2fc/3r3 _ 



1 - 



2M (i, r*) 



(111) 



where 



M(i, r*) = (Mo + AT*) - 



(go + g*r 

2r* 



(112) 



This differential equation governs the dynamics of the 
system in question, with different charge species and 
strengths giving different evolutions for the shell. It is in- 
teresting to note that for r* << 1 the exponential terms 
become extremely small (providing k ^ 0) and one finds 
dr+ldt w 0, implying r* 6 R. Hence, any inf ailing shell 
in this coordinate system will take infinitely long to reach 
r = 0. 

For completeness we write down the line element here 



dS 2 = 



2M\ ., 
- I 1 \dt 2 

r 



exp ( ^ ) dr 2 + rW, 



condition H103H is not satisfied. 



(113) 



where M(t,r) is given by equation (j9"8|) . One can see 
that the only function left to be determined in the line 
element is r* which appears implicitly through the mass 
function and is the solution of equation (jllip . 

Either side of the infinitely thin shell, the spacetime is 
vacuum and must necessarily be given by the RN space- 
time. Indeed evaluating the Einstein tensor for the line 
element given in (| 1 1 3[) either side of the shell gives an 
energy momentum tensor in accordance with the RN 
spacetime. Hence, the line element given above with 
M = M — Q 2 /2r (where M is a constant representing 
either M or Mq + M* and Q is either Q or Q + Q*) 
describes the RN spacetime. One can also look at the 
asymptotics of the spacetime and find that the metric 
tends to the Minkowski metric as r — > oo. This is con- 
sistent with the earlier description that this is equiva- 
lent to the marginally bound case with charge as the 
metric tending to the Minkowski metric implies the ob- 
server at infinity has zero energy. It is interesting to note 
that this spacetime is a specific form of that given by 
[26|, in particular, taking equation (13) from [26] with 
c(r) = exp(2fc/3r 3 ). 
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A. Apparent Horizon 



the metric to read 



An important facet of the discussion of specific solu- 
tions of equation (jlll|) and associated forms of the line 
element (|113|) is the evolution of the apparent horizon. 
It is well established that there are up to two horizons 
in the RN spacetime, which we denote as r±. In [33] we 
derived the apparent horizon for any spacetime with line 
element given by equation (f3"2")l with R = r to be 



(114) 



Substituting equation (|98|) into the apparent horizon and 
solving for r± in the two regions we find 



r ± (t) = 2M[t,r ± (t)}, 



r ± (t)=M ±\/Mo 2 -Qo 2 , 



(115) 



for r± < r* and 



r±(t) = M a + ± yJ{M + M*) z - (Q + G*) , 

(116) 

for r± > r± . These are the well known forms of the appar- 
ent horizon [ijj and we see there is a jump discontinuity 
across the thin shell of matter. 

We note here that the RN spacetime can be broken 
into four distinct subcases which has relevance for the 
apparent horizon and we will also see has relevance for 
the collapsing shells. Consider for the moment the RN 
spacetime with mass and charge M and Q respectively. 
The four cases are as such 

1. When Q = the spacetime is simply Schwarzschild. 
In this case r_ = and r + = 2M . 



2. When < Q < M the interior horizon is < r_ 
M and the exterior horizon is M < r + < 2M . 



< 



Q = M is known as an extremal black hole and we 
note that the interior and exterior horizons coin- 
cide, i.e. r_ = r + = M. 

A model of a classical charged particle is often con- 
sidered with Q > M. In this case we notice that 
the discriminant in the apparent horizon becomes 
negative and hence neither apparent horizon exist. 
This is therefore a naked singularity as the point 
mass at r — is visible by an observer at infinity. 



a 2 - 



fl 2 



2(3 dr 



1 



l + E l + Edt 1 + E \dt 



dr 



(117) 



where dr 2 — —dS 2 and we have reverted back to the 
3+1 notation for brevity in the expressions. We are only 
concerned with the motion of the shell and hence only 
interested in the relation of the proper time and coordi- 
nate time along the shell. Therefore, we can substitute 
equation (| 1 1 [) into the above, invert and take the square 
root to find the simple relation between the coordinate 
and proper times evaluated along the shell 



dt 

7h 



i 



■■ exp 



k 

3rJ 



(118) 



We note that when k = the coordinate time and proper 
time along the shell are equivalent. 

The proper motion of the shell can now be calculated 
from equation (jllip by substituting the above and we 
find 



dr* 

~d7 



/ e -2fe/3r; 5 



2M(r,r t ) 



(119) 



One can gain much information about the motion of the 
shell by looking at its acceleration. This gives insight into 
the relativistic "forces" in the spacetime. Differentiating 
with respect to the proper time we find 



«* 



3 -2fc/3rJ 



m , (Go + e*) 



(120) 



where a*(r) := d 2 r i ,/dr 2 is the proper acceleration of the 
shell and we have defined 



M := M + M*, 



(121) 



for convenience. It is of interest to explore the regimes 
where the acceleration is positive and negative as this 
gives us insight into the nature of the motion of the shell. 
However, the expressions for the velocity and acceleration 
of the shell are difficult to analyse in general, so we re- 
vert to specific examples to highlight different physical 
aspects. 



B. Shell proper motion 

One can study the motion of the shell in a coordinate 
invariant way by looking at the proper motion of the 
charged shell. The proper time is found by looking at 
the geodesic equation coming from the metric. Spheri- 
cal symmetry implies the motion of the shell is indepen- 
dent of the angular coordinates and one can rearrange 



C. Shell Examples 

1. Zero Charge 

As a first example we reduce the equations to those of 
non-charged dust. It is interesting to see that there is 
still a difference between a non-charged shell of matter 
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falling onto a Schwarzschild black hole as compared to an 
RN black hole. This is exhibited in the equation govern- 
ing the evolution of the shell (jllip as the Schwarzschild 
black hole has Qo = whereas it is nonzero for the RN 
black hole. Despite the non-charged shell not interacting 
through EM forces in the RN case, the geometry of an RN 
black hole is significantly different to the Schwarzschild 
black hole, which is exhibited through the different mo- 
tion of the shell. 

For this example we look at the evolution of the non- 
charged shell simply falling onto a Schwarzschild black 
hole, i.e. Qo = Q* = 0. This implies from its defi- 
nition that the constant k vanishes, which implies the 
coordinate time and the proper time are equivalent. Fur- 
thermore, equation (| 1 19[) greatly simplifies and can be 
integrated to give 



3/2 



r*(0) 



3/2 



(122) 



This equation is exactly that of a geodesic for a zero en- 
ergy particle and is the model elucidated by Adler et al. 
[331 ]. It is trivial to see from the above that as proper 
time evolves, r* moves closer to r — and will even- 
tually collide with the singularity, as expected. We can 
calculate the acceleration of the shell here such that we 
can compare with later results and we find 



M 



(123) 



As the mass is positive definite the acceleration of the 
shell in this case must always be negative. Therefore, 
as the only influences on the shell of particles are grav- 
itational, we discern that a negative acceleration corre- 
sponds to an attractive "force". 3 



2. Zero k 



Analytic solutions of (|119p can also be found when 
k = 0. We see this is the case when either 



Q* = and Qo £ 0, 



2Q =-Q* 



(124) 



(125) 



In fact, the dynamics of the entire spacetime are indis- 
tinguishable for these two charge distributions. Thus an 
uncharged shell of matter falling onto an RN black hole is 
equivalent to the charges being given by equation (|125|l . 
For the remainder of this section we shall discuss the 



Qi, = case, however we note this is equally applicable 
to the 2Qo = — Q* case. 

Putting k — into the line element we find it becomes 



dS 2 = - 1- 



2M 



dt 2 



2M 



dtdr + dr 2 + r 2 dn 2 



(126) 



where M is still given by (JM|) with = 0. This 
line element is the RN spacetime in Painleve-Gullstrand 
type coordinates both interior and exterior to the 
shell. While the Painleve-Gullstrand coordinates for the 
Schwarzschild spacetime cover the regions from r = to 
infinity, the equivalent version for the RN spacetime do 
not. One can see this by looking at the gt r component of 
the line element, which reads 



9ti 



I2M 
r 

We see that this coefficient becomes complex for 



Q2 
-2 1 



r <r c 



2M ' 



(127) 



(128) 



where we have defined the critical radius of the space- 
time, r c . That is, the coordinates are invalid for r < r c 
and therefore they have a range given by r e [r c , oo). 
One can compare this to the position of the horizon in 
the RN spacetime. In particular taking two extreme lim- 
its of the RN black hole, namely Q = and Q = M. 
Firstly, we find when Q = the critical radius also van- 
ishes and hence is not an issue. At the other limit to 
consider is that of an extremal black hole where Q = M 
one finds r± = M > M/2 = r c . Hence we deduce for 
all RN black holes such that < Q < M the critical 
radius is inside the apparent horizon. For black holes 
with Q > M the apparent horizon does not exist and the 
critical radius is therefore visible to an external observer. 

Possibly one of the most interesting things about 
this line element is the potential link with microscopic 
physics. The critical radius defined in (1128|) is actually 
the Compton radius which describes the "electron ra- 
dius" (see for example [1^]). Furthermore, we know that 
the apparent horizon, r + , for an electron is necessarily 
inside this radius, r + < r c . Hence, this line element is a 
way of classically describing the gravitational and electro- 
magnetic field of an electron without having to consider 
the complicated interior region or by modelling it as a 
point-mass. 

Reverting back to the collapsing shell within this space- 
time, the evolution of the shell becomes 



dr*_ 
dr 



I2M 



(129) 



3 This is actually due to the curvature of the spacetime, however 
we loosely use the term "force" as we will be comparing the 
effects due to gravitation with that of the Coulomb force. 



It is interesting that the sign of Qo does not matter to the 
evolution of this system. That is, changing the species 
of Qo from negative charge to positive charge does not 
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effect the system provided the sign of <2* also changes 
such that equation (|125p is satisfied. It can also be shown 
from equation (|129p that dr^/dr vanishes when 



r*(t) 



(130) 



Furthermore, when r* < r c the term in the square root 
sign of equation (|129|) becomes negative and hence there 
are no real solutions of this equation. 
Integrating equation (|129p gives 



C = 



-1 
3M 2 



2Mr*-Ql(Q 2 + Mr ic ), 



(131) 



where C is an arbitrary constant of integration. One can 
see when Qo = this equation reduces to equation (|122|) 
describing the non-charged shell as expected. By putting 
r* = r c into the above solution, one can show that the 
critical proper time, r c , defined to be the proper time at 
which the shell crosses the critical radius is given by 



1 



3M' 



Ql[Ql 



Mr A 



(132) 



which is positive definite. 

We can highlight specific features of the collapsing shell 
by looking at its acceleration throughout the spacetime 
as this is proportional to the "forces" acting on the shell. 
Evaluating equation (|120p with — wc find 



1 



Therefore we find 



a* < 



M 



M 



(133) 



(134) 



In direct comparison with the uncharged case (section 
I VII C ip . a negative acceleration corresponds to an at- 
tractive force on the shell. However, we see that for 
r c < r * < Qq/-M the acceleration of the shell is posi- 
tive and hence the force on the shell is repulsive. The 
value r a is therefore the point at which the acceleration 
changes sign. As <2* = in this case, the infalling shell 
is simply dust and hence there are no internal forces be- 
tween the particles in the shell. Furthermore, the shell is 
uncharged implying there is no Coulomb force from the 
charged point mass at r = acting on the shell. We saw 
that this repulsive term was not present in the case of an 
uncharged shell falling onto a Schwarzschild black hole. 
Therefore this repulsion is simply a manifestation of the 
geometry of the RN spacetime associated with the pres- 
ence of the extra term in the energy momentum tensor 
associated with the charge. 

We saw that the acceleration of the shell is positive for 
a finite range of the radius, however the position of this 
range as compared with the apparent horizon is inter- 
esting. Again we divide the spacetime into four distinct 
subcases; 



1. When Qo = we find r a = and hence the force 
is always attractive (this is the uncharged case dis- 
cussed in section I VII C ip . 

2. When Qo < Ad we can show that r_ < r a < r + . 
That is, the position at which the acceleration 
changes sign is hidden behind the apparent hori- 
zon of the RN spacetime. Therefore this effect is 
not observable from an external observer. 

3. When Qo = M. the two apparent horizons and the 
position of the change in acceleration also align, 
that is r± = r a = M. . 

4. The interesting case is again where Qo > M. im- 
plying the apparent horizons do not exist. In this 
case we find r a > M. , implying this effect is visible 
by an external observer. 



3. Exterior Schwarzschild 



As a next example we look at the case where 
Qo = -Q*. 



(135) 



This example is interesting as the mass function for the 
exterior region, r > r±, is simply M = Mo+M+, implying 
this exterior region is Schwarzschild. The charges in this 
case are equal but of opposite sign, which are therefore 
known to possess an attractive Coulomb force. We find 
from equation (| 104|) that k < and the evolution of the 
shell is given by 



dr* 
~d7 



I exp 
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-1 + ^1. (136) 



Furthermore, equation (|120p can be shown to be nega- 
tive definite implying the acceleration is always negative, 
hence the force is always attractive. 

Of interest in this model is that an observer external 
to the collapsing shell, i.e. sitting at some r > is 
in a region of spacetime equivalent to the Schwarzschild 
spacetime. However the observer can view the motion of 
the shell and infer the existence of the charges as the mo- 
tion of the shell is significantly different to the uncharged 
case. 



VIII. STATIC EQUATIONS 

Static charged fluid spheres have been studied in great 
detail, esp ecially perfect fluid spheres in diagonal coordi- 
nates [341 ] . The method we have utilized herein can also 
find static solutions by finding a timclikc Killing vector 
that is hypersurface orthogonal which implies the space- 
time is static and using this vector to reduce the 
general equations. Such a Killing vector is given by [l| 



(137) 
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Furthermore, solving Killing's equations with this vector 
field gives equations which further constrain the metric 
coefficients such that the spacetime is static. These equa- 
tions are given by 



=0. 



C n R - 

2{l + E) ' 
C n a =0 



-1B0 
a dr ' 



(138) 
(139) 
(140) 



Equations (|138|) and (|139[) together with equation ([35] 
imply that a static, spherically symmetric Einstein- 
Maxwell spacetime necessarily has vanishing heat flux 



3=0. 



(141) 



Putting the above conditions through the remaining non- 
trivial fluid and EM equations, ([36]), JSSJ), ([39]), ([54]) 
and ([53]) . we can derive a generalization of the Tolman- 
Oppenheimer-Volkoff (TOV) equations for hydrostatic 
equilibrium which include anisotropic stresses and elec- 
tric fields 



1 da 

a dr 



M + 4tt (P f - 2Tl F ) R 3 - Q 2 /2R OR 



R 2 (1 - 2M/R) 



dr 



(142) 



-1 



i F + P F - 2U F 



d_ 

dr 



P 1 



2W 



)_dR 
R~dr~ 



6U F - 



Qe 



R^/l - 2M/R / 
where the definition for the mass function is now 



dM 



dR Q 2 dR 



dr 4lTP R dr ' 2R 2 dr 
and the energy function is given by 

2 



(143) 



(144) 



(145) 



In the static limit, the electric field permeating the space- 
time is given by 



£ 



Qy/1- 2M/R dR 



R 2 



dr 



(146) 



It is trivial to show if we let R = r, then from (|138[) 
the shift necessarily vanishes, (3 = 0, and all dependancc 
on the temporal coordinate falls out of the equations, 
reducing to the familiar meaning of a static spacetime. 
Furthermore, reducing the energy momentum tensor of 
the matter to the perfect fluid case (II F = 0) results in 
the equations reducing to that of [36[. 



A. Exterior static region 

Analogously to the non-static cases, we can utilize the 
form of the metric for the interior regions of the spacetime 



to dictate the coordinates for the exterior, vacuum region. 
Thus, we again let all the energy-momentum terms van- 
ish at some finite radius, and find the charge and mass 
respectively become Q = Qo and M = Mo — Qo 2 /2R. 
Putting this into the equation of hydrostatic equilibrium 
(|142p . we find the lapse function can be integrated to give 



a= \ I - 



2Mn 



R 



Qo 2 
R 2 



=£r> ( 147 ) 



where coordinate freedom has been used to scale an arbi- 
trary function of the temporal coordinate to unity. There 
still exists coordinate freedom in the exterior which is de- 
termined from a more precise stipulation of the interior. 
For example, if we find a solution of the interior region 
that has R — r, this implies through equation (|138[) that 
(3 = 0. Equation (|145|) then implies E is known, and the 
metric is simply the RN metric. Therefore, a static inte- 
rior charged fluid, with R = r, naturally matches onto an 
RN spacetime in RN coordinates, such that both regions 
of the spacetime are expressed using the same coordi- 
nates. 



IX. CONCLUSION 

In this article, we sacrifice mathematical aesthetics 
for physically interpretable results, and show that many 
known solutions of the Einstein-Maxwell field equations 
that are often treated separately can be derived from the 
same ansatz. Furthermore, we showed that the matching 
of various regions of the spacetime can be done in single 
coordinate patches as they are all manifestations of the 
same underlying equations. We discussed the construc- 
tion of a qualitative example whereby an initial spacelike 
hypersurface was established with a general charged fluid 
interior and vacuum RN spacetime exterior. This system 
is evolved using the derived equations and a region of 
charged Vaidya emerges in a region between the charged 
fluid and the vacuum. Therefore, on any t equals con- 
stant slice subsequent to the initial hypersurface there 
exists a region of charged fluid ball extending to some 
finite radius with charged Vaidya beyond this and vac- 
uum beyond this extending out to spacelike infinity. This 
model provides causality to the charged Vaidya solution. 

The construction of these solutions using single coordi- 
nate patches allowed for the discussion of the collapse of 
thin shells onto an RN black hole. A surprising result of 
this was that a shell of charged matter with the opposite 
charge to the black hole passes through a region where 
the "force" on the shell was repulsive. This was shown 
to be attributed to the geometry of the RN spacetime. 

This work opens avenues for further research into many 
currently unanswered questions. In particular the study 
of the formation and evolution of shell crossing and shell 
focussing singularities. While this has been extensively 
researched for simple fluids models in spherical sym- 
metry, very little is understood about the inclusion of 
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charged fluids. This work will require numerical anal- 
ysis as one can quickly see that the equations derived 
herein are highly coupled, implying closed form solutions 
are generally not available. However, the formulation 
used herein implies this numerical treatment is achiev- 
able as the equations are all expressed as an initial value 
problem. In particular, in the case of thin charged shells 
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